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The Dusty Aftermath of the V1309 Sco Binary Merger 
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ABSTRACT 

We present mid-IR photometry and spectra of the merged binary V1309 Sco taken 
between 18 and 23 months after outburst. Strong mid-IR emission and a sohd state 
absorption feature indicate the presence of a significant amount of dust in the circum- 
steUar environment. The absence of detectable mid-IR emission before the outburst 
suggests this dust was produced in the eruptive merger event. 

Model fits to the solid state absorption feature constrain the constituent species 
and column density of the dust around V1309 Sco. We find the absorption feature can 
be reproduced by large (3 |i.m) amorphous pyroxene grains at a temperature of 800 K. 
This grain size, if confirmed with longer wavelength spectroscopy and modelling, would 
be suggestive of dust processing in the circumstellar environment. The data in hand 
do not allow us to discriminate between disk or shell configurations for the dusty 
material. 
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1 INTRODUCTION 

V1309 Sco is the firs t documented case of a binary merger 
l|Tvlenda et al.ll201lf ). Disco vered as Nova Sco 2 008 when it 
erupted in September 2008 (|Nakano et al.l [20081 ). its subse- 
quent evolution marked it as a new member of the 'red nova' 
class of stellar outbursts, with a spectrum rapidly evolv- 
ing from F to M spectral types (Rudy et al. 2008a, 2008b) 
and line profiles transi tioning from absorption to emis- 
sion llMason e t al?'2010l). The red novae were proposed by 
ISoker fc Tvle nda (20031) to be binary merge rs, and include 
the e qually enigmatic objects V838 Mon (|Tvlenda et al.l 
|2005^ and V4332 S g r (IKamiriski et al.ll2010l ). 

[Tvlenda et al] (|201ll ) presented the crucial evidence 
that revealed the nature of V1309 Sco, and possibly the 
class to which it belongs. They showed that the photomet- 
ric variability found in its OGLE /-band light curve clearly 
marks V1309 Sco's progenitor as an early K giant contact 
binary with an orbital period of ~ 1.4 days. This orbital 
period decreased in the 6 years prior to eruption, simultane- 
ous with a slow increase of the system brightness, and they 



found no evidence of periodicity in V1309 Sco's light curve 
after its eruption. They calculated that the energy released 
in the outburst was consistent with the dissipation of or- 
bital energy expected for a merger between the progenitor's 
components. 

Assuming V1309 Sco's outburst was indeed the result 
of a stellar merger, questions remain as to the nature of the 
remnant. In many binary models a phase of deep contact 
is reached as the orbit shrinks, resulting in mass and angu- 
lar momentum loss from the system. This phase may lead 
to the formation of an excretion disk, whic h could remain 
around the coalesced star after the outburst (|Webbinklll976l : 
llben fc Liviolll993l : [Shu et al.lll979l ). Mid-IR measurements 
of the remnant are an ideal way to probe dust production 
during the merger process. 

In this Letter we show that V1309 Sco has become dom- 
inated by mid-IR emission since eruption. This can be at- 
tributed to dusty circumstellar material that likely formed 
during the merger, and we comment on the dust species that 
may constitute this matter. We also present evidence of dust 
processing. 
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Table 1. Mid infrared photometry of V1309 Sco 
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Sil 
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18.3 
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2 OBSERVATIONS 

The Optical Grav itational Lensing Experiment (OGLE; 
lUdalski et aLlll992[ l provides long-term optical light curves 
of many sources in the Galactic centre and the Magellanic 
clouds. V1309 Sco was monitored during the OGLE-III and 
OGLE-IV phases of the project, resulting in thousands of 
regular photometric measurements spanning more than ten 
years. The OGLE III + IV 7-band hght curve of V1309 Sco is 
shown in Fig. [T] Here we include new photometry from the 
2011 and 2012 observing seaso ns of the OGLE-IV project 
that was not presented in Tvlenda et al.l (|201lh . 

The Wide-fie ld Infrared Survey Explorer (WISE; 
I Wright et ai]l2010[ ) satellite surveyed the entire sky in four 
mid-IR photometric bands centred on 3.4, 4.6, 12, and 
22 nm. V1309 Sco was observed on 2010 March 19, with 
detections in all four bands. We downloaded the magni- 
tudes from the NASA/IPAC Infrared Science Archive and 
converted th e m to fluxes in Jansky using the formula in 
IWright et all l|2010l ). and applied colour corrections appro- 
priate for a 400 K blackbody (see Section [3.2|) . We do not 
report the more recent WISE 3-band Cryo Release photom- 
etry of V1309 Sco as its reliability is unclear at this time. 

V1309 Sco was observed by the Thermal-Region Cam- 
era Spectrograph (TReCS) at Gemini South under program 
GS-2010B-C-7 on 2010 August 3, as part of a project to char- 
acterise dust forming novae. We downloaded the raw N-band 
spectrum and narrow band images from the Gemini Sci- 
ence Archive. We extracted photometry of V1309 Sco using 
standard IRAF procedures and converted the magnitudes to 
fluxes in Jansky. We assume 10 per cent error bars on these 
measurements. We reduced the mid-IR spectrum using the 
Gemini IRAF mid-IR spectral reduction pipeline, msreduce, 
and flux calibrated by shifting it vertically to match the flux 
levels of the narrow-band photometry. 

The WISE and TReCS mid-IR photometry is shown 
in Table [1] The photometry and the TReCS spectrum are 
plotted in Fig[l 

3 RESULTS AND DISCUSSION 
3.1 Optical Light Curve 

The features that lTvlenda et all (j201ll ) noted in the OGLE 
7-band light curve, namely the slow brightening of the pro- 
genitor, the fading prior to eruption, the smooth outburst 
beginning ~ 6 months before its discovery as Nova Sco 2008, 
and the decline to magnitudes fainter and with smaller scat- 
ter than the progenitor, can be clearly seen in Fig. [T] With 
the addition of more recent photometry, it is clear that 
V1309 Sco for the most part continues to decline in bright- 



ness, opposite behav iour to that shown by red novae V 838 
Mon and V4332 Sgr (|Tvlendall2005l : iTvlenda et al.ll2005l ). It 
remains to be seen whether the brightening starting at JD 
~ 2456100 days is periodic or the beginning of a reversal. 
There also appears to be a slight increase in scatter over the 
two most recent observing seasons. 

Using IRAF 's Phase Dispersion Minimisation tool 
(|Stellingwerj|l978l ). we find no evidence of periodic variabil- 
ity in the new photometry. This lack of periodicity indicates 
that the variability of the remnant is unlikely to be pulsa- 
tion or binary-induced variations. It is unclear at this stage 
whether the post-outburst scatter and faintness is an intrin- 
sic result of the merger. We discuss an alternative scenario 
in Section [3.31 



3.2 Mid-Infrared Measurements 

We searched for previous mid-IR measurements of V1309 
Sco using its WISE coordinates of 17:57:32.94, -30:43:09.84 
which are accurate to 0'.'043 in RA and 0'.'044 in Dec. No 
detection was found at this position in the IRAS or AKARI 
catalogues. At 18 |xm, AKARI is 100 per cent complete down 
to 0.3 Jy (Ishihara et al. 2010), so at that epoch (2006) 
V1309 Sco was fainter than 0.3 Jy at 18 iJ.m, a difference of 
more than 9 times the 2010 TReCS Qa-band measurement 
(Table [l|. We note that AKARI detected a bright source 
(2.7 Jy) in this band 2 arcmin away from the position of 
V1309 Sco, further suggesting AKARI should have been ca- 
pable of detecting V1309 Sco if it was as bright as it was in 
2010. At 12 \im, IRAS is estimated to be 95 per cent com- 
plete down to 0.46 JjQ- From its non-detection at the epoch 
of IRAS (1983), V1309 Sco was likely fainter than 0.46 Jy at 
12 ^m, a difference of more than 6 times the 2010 TReCS 
Si5 flux. 

Although V1309 Sco was undetected in the mid-IR be- 
fore its eruption, as of 2010 it had reached flux densities of 
over 3 Jy at some mid-IR wavelengths (as seen in Table \T\ 
and Fig. [2]), despite the optical flux having declined to pre- 
outburst values. This suggests that the merger between the 
progenitor's binary components produced not only a huge 
optical flare but also a significant amount of dust. The pre- 
burst fading of the syste m (visible at JD ^ 24 54500 days in 
Fig. [T} was proposed bv I Tvlenda et af] l|201ll) to signal the 
formation of a dusty excretion disk. The recent high mid- 
IR flux certainly supports the hypothesis that the binary 
merger triggered signiflcant dust formation, and that this 
dust has been retained in the circumstellar environment of 
the system. 

To learn more about the dust surrounding V1309 Sco 
we made a simple model flt to the mid-IR measurements. 
The model takes the form of 



E 



{aiAiB„,i) 



(1) 



where D is the dust model, / is a constant scaling the total 
received flux, Bi^^cont is a blackbody representing the dust 
continuum emission, a; is a constant scaling the amount of 
the ith absorbing species, Ai are the absorption coefficients 
for the ith absorbing species, and B^^i is a blackbody at the 



^ IRAS Explanatory Supplement, 

|http://irsa.lpac.caltech.'edu/IR ASdocs/exp.sup/| 



© 2013 RAS, MNRAS 000, [HS] 



The Dusty Aftermath of the VI 309 Sco Binary Merger 3 




2500 3000 



3500 4000 4500 

JD - 2450000 (days) 



5000 5500 6000 6500 



Figure 1. OGLE /-band light curve for V1309 Sco. The peaic luminosity is not covered as the object was overexposed in the OGLE images. 
Gaps denote times when the object was too close to the Sun to observe. The dashed line shows the epoch of the WISE observations, 
and the dot-dashed line shows the epoch of the TReCS observations. The dotted line shows the approximate epoch of the AKARI 
non-detection; see Section 13.21 



temperature of ith absorbing species. Due to the sparseness 
of the data we allow only a single dust continuum, and one 
absorbing component. 

W e used absorption coefficients calculated bv lMin et al.l 
(2007'), using both Gaussian Random Field (GRF) and Dis- 
tribution of Hollow Spheres (DHS) grain shape distributions. 
Our dust species grid included amorphous olivine and py- 
roxene, amorphous Mg-rich olivine and pyroxene, crystalline 
olivine and pyroxene, and silica. For each species various 
grain sizes between 0.1 and 4.0 |J.m were considered. Our 
temperature grid covered T — 200 K to 1500 K in 50 K steps. 
We determined the best fit to the data by minimising ■ 

We find the dust continuum is best fit with a tem- 
perature of 400 K. Of the features in V1309 Sco's mid-IR 
spectrum, most notable is the broad absorption band at 
9.7 [im, usually attributed to amorphous silicates. We find 
the best fit to the absorption features is amorphous pyrox- 
ene MgFeSi206 grains 3.0 \im in size, at 800 K, with a DHS 
shape distribution. In general we found that was more 
sensitive to grain size than to species, temperature, or shape, 
with larger grain sizes (3.0 and 4.0 vim) corresponding to 
lower values, reproducing the broad 9.7 \J.m feature with 
greater accuracy. The best fit is shown in Fig [2] 



It is clear that even our best fit is somewhat lacking. 
This poorness of fit is likely due to the simplistic nature of 
our model, however a more complete SED is required before 
detailed modelling is undertaken. Although we quote the 
best fit species, temperature, and shape distribution as likely 
properties of V1309 Sco's circumstellar dust, we note here 
that it is difficult to properly constrain these with the data 
in hand. However we are confident that the broad 9.7 [im 
feature indicates large grains. 

Because V1309 Sco's mid-IR SED shows dust absorp- 
tion instead of emission features, indicating an optically 
thick circumstellar environment, we are able to derive only 
the column density of the dusty circumstellar matter, not 
the total dust mass. Column density may be calculated by 
dividing the optical depth by the extinction cross section 
of the absorbing grains. The optical depth at a particu- 
lar wavelength can be calculated from the observed fiux 
and the continuum model fit; the extinction cross section is 
c"A = QxT^o? , where Q\ is the extinction coefficient of that 
particular grain type and size at a particular wavelength, 
and a is the grain radius. We used extinction coefficients for 
'astrono mical silicate' a t A = 9.7vim and a = 1.5[i.m calcu- 
lated bv lLaor fc Drain3 (|l993l ), as extinction coefficients for 
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Figure 2. Mid-IR SED of V1309 Sco. The TReCS spectrum is shown by the black hne, the WISE photometry by the green squares and 
the TReCS photometry by the blue diamonds. The red line represents our simple dust model fit. The dashed purple line is a blackbody 
fit to the continuum at a temperature of 400 K. The WISE error bars are comparable to the point size. 



our modelled dust grains were unavailable. We find that in 
the 9.7 [im feature, the column density of the 3 \im amor- 
phous pyroxene grains is 4 .3 x 1 0^ cm~^. 

As iBaird fc Cardellil (|l985l ) note, column density and 
total dust mass may not necessarily correlate, as an extended 
dust cloud may have a large mid-IR excess (indicating high 
mass) but low column density due to its large volume. There- 
fore without an estimate of the radial extent of the dust, we 
are unable to estimate its mass. We stress that the column 
density derived here is a lower limit for the total dust col- 
umn density of V1309 Sco's circumstellar environment, as it 
includes only the 3 \im absorbing pyroxene dust grains, and 
neglects any particles that may contribute to the dust con- 
tinuum emission. The details of V1309 Sco's circumstellar 
environment are shown in Tabled] 

Two particular results should be noted from our dust 
fitting. The first is that the warm temperature of the dust 
suggests it formed recently, consistent with the absence of a 
previous mid-IR detection. The second is that the large size 
of the grains producing the absorption features suggests the 
dust has been processed in the circumstellar environment: 
dust formed in the atmospheres of AGB stars and returned 
to the ISM is typically sub-micron in size, and larger grains 
require processing in a stable environment. It is difficult to 



Table 2. The Circumstellar Environment of V1309 
Sco 



Parameter Value 

^cont 

400 ±25 K 

Tfcat 800 ±25 K 

Dust species Amorphous pyroxene MgFeSi2 06 

Grain size '^3.0 \ira 

Column density" (4.3 ± 0.43) X cm'^ 

Optical depth" 0.71 ± 0.07 



" Calculated at the 9.7|J.m feature 

say at present whether V1309 Sco's dusty shroud is a disk 
or a shell, but the processing of dust grains to sizes of a few 
microns certainly implies a quick evolution, if we assume the 
dust was formed as part of the recent merger event. 

3.3 The Nature of V1309 Sco 

From the work of lTvlenda et al.l (|201lh and the results pre- 
sented here, it seems that V1309 Sco was a contact binary 
that produced a nova-like outburst as the two stars merged, 
and is now surrounded by warm optically thick dusty matter 
showing strong absorption from processed amorphous pyrox- 
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ene particles. The lack of periodicity in the post-outburst 
lightcurve supports the hypothesis that V1309 Sco is now a 
single star. It seems likely that the dusty envelope formed 
during the merger. 

Many models of binary interaction involve mass loss and 
formation of a disk, such as the circumstellar excretion disk 
proposed bv lWebbinkl l 19761 ). Because V1309 Sco's progen- 
itor was a photometrically varying contact binary, we are 
viewing the system ne ar its equa torial plane, but not neces- 
sarily exactly edge on l|Nicholls fc Wood 2012 1 . It is possible 
that the large amount of dust we detect resides in a disk seen 
slightly tilted from edge-on, meaning although it produces 
the large mid-IR excess we observe it does not completely 
enshroud the system, extinguishing only part of the optical 
flux. This could account for the recent faint /-band magni- 
tudes of V1309 Sco as seen in Fig. [Jl and the pre-outburst 
fading. The post-outburst irregular small variations seen in 
the light curve may also be caused by a slightly tilted disk, 
as clumps of dust pass in front of the line of sight. How- 
ever, with our current data we are unable to discern whether 
V1309 Sco's dust is in a disk or another configuration. 

A similar mid-IR excess was recen tly observed in y 4332 
Sgr, a red nova that erupted in 1994 (jSaneriee et al.| [2007!: 
iKamihski et al.ll2010l 'l. It too shows a large 9.7 ^.m absorp- 
tion feature indicating an optically thick ox ygen rich dusty 
enviro nment, and evidence of alumi na. Both Baneriee et al.l 
l|2007ll (and references therein) and IKamihski et al.l "( 2010l l 
proposed that V4332 Sgr's dust occupies a stable disk that 
is probably edge on to our line of sight and hence causing 
the observed high IR-to-optical luminosity ratio. No previ- 
ous measurement of the column density of V4332 Sgr's dust 
has been made, but we can make a rough estimate if we as- 
sume the grains are similar to what we find for V1309 Sco. 
In this way, we estimate that the column density of the ab- 
sorption feature seen toward V4332 Sgr is ~ 5 x 10^ cm"'^, 
comparable to the value we calculate for V1309 Sco. More 
detailed dust modelling would constrain this value further. 

V1309 Sco, as the first confirmed binary merger, pro- 
vides an ideal laboratory to test merger models. Proposed 
merger remnants in the literature include blue straggl ers, 
rapid ly rotating giants like FK Com (Gazcas & Stc piehl 
l2008h. and giant s accreting from a circumstellar disk 
l|Melis et aLll2009t ). It remains to be seen whether V1309 
Sco shows properties indicative of these types of objects. 

V1309 Sco is a most fascinating object and worthy of 
further study. Continued photometric monitoring is impor- 
tant to track its luminosity and variability evolution. Optical 
spectra are necessary to provide measurements of its rota- 
tional velocities. Mid-IR spectra at longer wavelengths are 
essential to determine the exact composition of its dusty cir- 
cumstellar environment, and mid-IR monitoring will reveal 
any time evolution of this dust. Whether the dust surround- 
ing V1309 Sco resides in a disk is, we feel, one of the most 
pressing questions. 
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